Dynamics of dispersive single qubit read-out in circuit quantum electrodynamics 
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The quantum state of a superconducting qubit non-resonantly coupled to a transmission line 
resonator can be determined by measuring the quadrature amplitudes of an electromagnetic field 
transmitted through the resonator. We present experiments in which we analyze in detail the dy- 
namics of the transmitted field as a function of the measurement frequency for both weak continuous 
and pulsed measurements. We find excellent agreement between our data and calculations based 
on a set of Bloch-type differential equations for the cavity field derived from the dispersive Jaynes- 
Cummings Hamiltonian including dissipation. We show that the measured system response can be 
used to construct a measurement operator from which the qubit population can be inferred accu- 
rately. Such a measurement operator can be used in tomographic methods to reconstruct single and 
multi qubit states in ensemble averaged measurements. 

PACS numbers: 42.50.Ct, 42.50.Pq, 78.20.Bh, 85.25.Am 
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I. INTRODUCTION 

Among several stringent requirements like scalability 
or precise coherent control, high-fidelity read-out of the 
qubit state is an important aspect of all experimen- 
tal efforts in quantum information science [1]. For su- 
perconducting qubit s l2|l a number of read-out strate- 
gies jH, 0, [H, 0, 0, H, IU Ej specific to various implemen- 
tations have been pursued. Early charge qubit read-outs 
implemented with single electron transistors (SET) 
were limited by strong current noise back-action from the 
measurement device on the qubit. Similar limitations 
applied to flux qubit read-out, measuring the switch- 
ing current of a nearby superconducting quantum in- 
terference device (SQUID) [H, [l3[. One strategy to 
achieve high read-out fidelities is to perform a quantum 
non-demolition (QND) measurement whichpreserves the 
eigenstates of the system Hamiltonian [14[. Repeated 
measurements yield identical results and consequently 
an improved signal-to-noise ratio. In the circuit quan- 
tum electrodynamics (QED) architecture, where a su- 
perconducting qubit is strongly coupled to a transmis- 
sion line resonator [TBI, [l6[ , the qubit can both be con- 
trolled and read-out via the cavity using microwave sig- 
nals. The read-out can be accomplished by detecting the 
dispersive qubit-state dependent shift of the resonator 
frequency [13] • In the dispersive limit, where the qubit 
transition frequency is far detuned from the resonator 
frequency, and for small photon numbers the measure- 
ment of the transmitted resonator field forms a QND 
measurement [HI, EE EE dEH. Similarly, QND mea- 
surements have been employed to demonstrate explic- 
itly the repeatability of this type of measurement for a 
flux qubit dispersively coupled to a nonlinear oscillator 
[20l I21L |22| . In a related experiment the phase of a tank 
circuit coupled to a flux qubit was monitored, demon- 
strating resonant tunneling [23j]. Note also, that in the 
circuit QED architecture a QND measurement has been 
proposed to g enerate and detect multi-qubit entangled 
states (H [lalli. 



Here we analyze the time-dependent response of the 
quadrature amplitudes of an electromagnetic field trans- 
mitted through a resonator to changes in the qubit state 
under dispersive interaction. We derive and discuss a 
set of Bloch-type equations describing accurately the dy- 
namics of the qubit and the resonator for continuous mea- 
surements. We also extend the analysis to pulsed read- 
out which avoids measurement-induced dephasing during 
qubit manipulation and allows for stronger measurement. 
We analyze experimental data for different measurement 
frequencies and find excellent agreement with theory us- 
ing a single set of independently measured parameters. 
This approach has already been successfully used in ex- 
periments pj], H3, HE but not yet discussed in litera- 
ture. 

Within this framework, we demonstrate our ability 
to infer the state of the qubit embedded in the cavity 
from a measurement of both field quadratures transmit- 
ted through the resonator. The construction of the corre- 
sponding projective measurement operator based on the 
state-dependent resonator response is outlined. For Rabi 
oscillation measurements, we discuss the extraction of the 
qubit state population and compare to numerical simu- 
lations. 



II. THE PHYSICAL SYSTEM AND ITS MODEL 

We consider a transmon-type qubit [lO, HH embedded 
in a transmission line resonator, as illustrated schemat- 
ically in Fig. [TJ The qubit is coupled to the resonator 
through the effective capacitance Cg, le ading to a qubit- 
resonator coupling of strength g |l6u32| . The qubit tran- 
sition frequency cj a is tunable by an externally applied 
flux [30[. The resonator with resonance frequency u r 
determined by its geometric and dielectric properties [33| 
is modeled as an LC circuit. 

When the qubit /resonator detuning A ar = cj a — uo T is 
much larger than the coupling strength g, this system is 
described by the dispersive approximation of the Jaynes- 
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FIG. 1: (color online) Circuit diagram of the experimen- 
tal setup. A harmonic oscillator modeled as an LC circuit 
with resonance frequency oo r is coupled to a transmon-type 
qubit [30|] through the coupling capacitance C g . The qubit 
transition frequency co a is controlled by an externally applied 
magnetic flux The qubit state is coherently manipulated 
by a pulsed microwave source at the frequency uj s . The res- 
onator is probed by a signal applied to the input capacitor d n 
at the frequency cj m . The transmitted signal is amplified and 
downconverted by mixing with a local oscillator at frequency 
cjlo and then digitized using an analog to digital converter 
(ADC). 
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is the dispersive coupling strength between the resonator 
and the transmon qubit approximated as a two-level sys- 
tem [30]. E c is the charging energy. The dispersive 
coupling leads to a qubit state-dependent shift of the 
resonator frequency, which we use to measure the qubit 
state. 

As illustrated in Fig. [TJ the qubit state is controlled 
by a coherent microwave field of amplitude Q(t) and fre- 
quency cj s applied directly to the qubit while the mea- 
surement tone with amplitude e m (t) and frequency co m is 
applied to the input port of the resonator. These exter- 
nally applied control fields are modeled by the Hamilto- 
nian 



H d ft (e m (£)a T e 



ft(£)<7 + e- ia;st + h.c.), (3) 



where we have taken e m (t) and Q(t) to be real for sim- 
plicity of presentation. 

The dynamics of the system in presence of dissipation 
and dephasing is described by a Lindblad-type master 
equation [34| 



p = -^[^,p] + ^[a]p+ 7l P[a_]p+^P[a z ]p^/:p, (4) 



where H = H disp + H d and V[A]p = ApA^ - A^Ap/2 - 
pA^ A/2. Here, 71 = 1/Ti is the qubit decay rate, 7^ the 
qubit pure dephasing rate and k, the photon decay rate. 
Since all experiments discussed in this paper are done at 
a small photon number n <C n ncr j t = |A5; r |/4g 2 and as 
71 exceeds the Purcell decay rate [35[ we neglect higher- 
order corrections to this dispersive master equation [36|, 
[33. 

To study the dynamics of the coupled qubit /resonator 
system, we derive Bloch-like equations of motions for the 
expectation value of the qubit operators (&{) (i = x, y,z) 
and the resonator field operators (a) and (aJa). However, 
the master equation (|4j) leads to an infinite set of coupled 
equations for these expectation values. For instance, the 
differential equation for (a) involves terms proportional 
to (a& z ), (a^aacJz) and (a<r x ), which in turn involve even 
higher order terms. We therefore truncate this infinite 
series by factoring higher order terms (a^a&i) ~ (a^ &)(&[) 
and {a^aa&i) w (aJa)(aa[), but keeping the terms (a&{) 
which ensures that the field contains information about 
the qubit state. This choice of factorization yields the 
correct average values for coherent and Fock states [381 ] 
and leads to a complete set of eight coupled differential 
equations 



d t (a) 
dt(crz) 

dt(a y ) 

dt(aa z ) 
dt(aa x ) 
dt(aay) 

d t {a) a) 



i A rm 

n<* y > - 



(11 
" V 2 

-i A rm 



(a) - ix(aa 2 ) - 
- 71 (1 + (&.)) , 
+ 2 X (Va) + 

+ 7*) (o"x), 

-2x(Va> + ^ 

+ 7*) (CTy) - 0(<X Z ), 

(d<r z ) - ix(a) + n(dd-y) 
(a z ) -71(d) - (71 



■<&>> 



<*y> 



(dx) 



(d«T z 



(5a) 
(5b) 

(5c) 

(5d) 
(5e) 



iA rm (d<j x ) - [A £ 
(d x ) - ( : 



-ie 



'7i 
2 



-iA rm (aa y ) + [A a 

-^mK) - (^y +70+ 2 J \ 6 
-Q(a(j z ), 

-2e m Im(a) — ft(a^a), 



+ 2 X ((a t a) + l)] (aa y ) 
(aa x ), (5f) 
1)] <aa x > 



70 



ft 

+ 2 

2 X ((d f d) 



(5g) 
(5h) 



which we refer to as Cavity-Bloch equations. Here, we 
have defined A as = co a — co s and A rm = uo T — uo m as 
the detuning of the control and measurement microwave 
fields from the qubit and cavity frequency, respectively. 
While these equations are apparently more complex than 
Eq. (j4j) , they can be analytically solved in some cases and 
are much faster to solve numerically. Note, that they do 
not include measurement-induced dephasing caused by 
photon shot-noise [15|, [39[ , because only the expectation 
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value of a) a is taken into account, and higher order mo- 
ments are omitted. This is of no consequence for the un- 
derstanding of the experiments presented here because 
of the small photon numbers present during the mea- 
surement. In experiments where measurement-induced 
dephasing is important, Eq. (|4]) has to be solved di- 
rectly [19f, or higher order terms have to be taken into 
account [361 ]. 



III. CONTINUOUS MEASUREMENT 
RESPONSE 

To experimentally determine the state of the qubit, 
we probe the dynamics of the resonator-qubit system 
by measuring the resonator transmission at different fre- 
quencies co m . A time resolved, phase sensitive measure- 
ment of the transmission quadrature amplitudes is re- 
alized by down converting the measurement signal at 
frequency co m in a mixer to an intermediate frequency 
A m LO = — ^lo = 25 MHz using a local oscilla- 
tor of frequency c^lo and phase 0, yielding one inde- 
pendent data point every 40 ns, see Fig dJ This inter- 
mediate frequency (IF) signal is digitized in an analog- 
to-digital converter (ADC). The whole experiment is 
repeated and averaged 650'000 times to enhance the 
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FIG. 2: (color online) a) In a continuous measurement, the 
cavity is driven with the qubit in its ground state, populating 
the resonator with n — 1 photons on average at the cavity 
resonance frequency. The qubit is then prepared in the ex- 
cited state with a 7r-pulse ending at t = 0. b) and c) averaged 
measurement response Q, / versus time t for a continuous 
weak measurement at the frequency u m = cj r — x- Solid lines 
show the predicted response from the Cavity-Bloch equations, 
Eq. (|5]). Time resolved data taken at different detunings A mr 
is shown in d) and e). The arrows indicate the detuning at 
which the data shown in b) and c) is taken. The colormap 
codes red for a positive amplitude and white for zero. 



signal-to-noise-ratio and then digitally down converted 
to DC, leading to a measurement time of 26 ms per 
point for a total of 6.5 seconds for 250 points in the 
trace. From this we obtain both the in-phase (/) and 
the quadrature (Q) components of the transmitted field 
Asin(A mL o£ + 0) = /sin A mL o^ + Qcos A mL o^ Using 
input-output theory [14( , these quadratures at the output 
of the resonator are related to the cavity Bloch equations 

by 



I(t) = ^Zhw.n Re(a(t)}, 
Q(t) = sJZhbj.K lm(a(t)), 



(6) 



where Z is the characteristic impedance of the transmis- 
sion line connected to the resonator. 

When arbitrary qubit rotations can be performed, it is 
sufficient to consider the measurement response for the 
qubit prepared in either its ground \g) or excited state 
|e) for a full characterization of the qubit state [4pj |. Fig- 
ure [2^) shows the pulse scheme used for the measure- 
ment. The time dependent quadrature amplitudes / and 
Q are measured at the cavity resonance frequency with 
the qubit in the ground state (c^ m = u Y — x). The res- 
onator is continuously driven at a measurement drive am- 
plitude of = k;/2, populating the resonator with n « 1 
photons on average in resonance. A 10 ns long tt pulse 
ending at time t = and resonant with the ac- Stark 
(4lj and Lamb shifted [42[ qubit transition frequency 
= [^a + 2x(( a ^ a ) + 1/2)] := ^s,res is then applied to 
the qubit, see Figs. [2)3) and c). Qubit relaxation during 
the 7r pulse limits the achievable \e) state population to 
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FIG. 3: (color online) Transmission spectrum of the resonator. 
The Q quadrature of the field is shown in a) and the / quadra- 
ture in b). The datapoints represent the instantaneous aver- 
aged response of the field before (red solid points), 180 ns 
(blue diamonds) and 740 ns after (green crosses) the 7r-pulse 
is applied. The underlying lines show the numerical simu- 
lations and the dotted line shows the expected response of 
the system for the qubit in the excited state |e) with infinite 
lifetime. 
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99% [43] , as obtained by solving the Cavity-Bloch equa- 
tions. This is within the statistical uncertainty of the 
detection. Furthermore, thermal excitations of the qubit 
are expected to be very low and are therefore neglected. 

The dependence of the quadrature components / and 
Q on the detuning A mr of the measurement frequency 
from the bare resonator frequency is plotted in Figs. [2]i) 
and e). For clarity, the quadratures are rotated in the 
/Q-plane for each measurement frequency cj m such that 
the Q quadrature is maximal in the steady-state (qubit 
in the ground state), resulting in Q = A and I = for 
t — > oo. As a result, before the 7r-pulse the / quadrature 
is always zero. 

The time and frequency dependence of the measure- 
ment signal is accurately described by the Cavity-Bloch 
equations with a single set of independently measured, 
non adjustable parameters as indicated by the solid lines 
in Figs. [2}d) and c). The cavity resonance frequency is 
determined as uu T /2n = 6.44252 ± 0.00002 GHz with a 
photon decay rate of ^/2tt = 1.69 ± 0.02 MHz. The 
qubit transition frequency is determined spectroscopi- 
cally as cj a /27r w 4.009 ± 0.001 GHz with a charging 
energy of E c /h = 232.5 ± 0.5 MHz [3l|. The transition 
frequency is adjusted using external magnetic flux. The 
qubit-cavity coupling g/2ir = 134 ± 1 MHz is extracted 
from a measurement of the vacuum-Rabi mode splitting 
at cj a = uo v [l6[ . 

The cavity pull x/2tt = -0.69 ± 0.02 MHz is de- 
termined spectroscopically. This is done by measuring 
the cavity resonance frequency leaving the qubit in the 
ground state and then measuring its frequency shift ap- 
plying a continuous coherent tone to the effective qubit 
transition frequency uj Syres . When the qubit transition 
is saturated (Q ^> 71), the resonator is shifted on aver- 
age by x- This value is in good agreement with the full 
transmon model taking into account higher levels [3(j 
(x/2tt = -0.71 MHz). 

In fitting the measurement response in Fig.[2j the qubit 
decay rate 7i/27r = 0.19 ± 0.01 MHz is used as an ad- 
justable parameter which is equal to an independent mea- 
surement of 71 within the statistical uncertainty. In prac- 
tice, the qubit decay rate is determined for one measured 
trace, and then kept fixed for all other traces. Note that, 
for short 7r-pulses, the dephasing rate 7^ has no measur- 
able influence on the solution of the equations. Addi- 
tionally, a single scaling factor is introduced to relate the 
quadrature voltages at the output of the resonator to the 
digitized voltages after amplification. 

To interpret the time and frequency dependence of the 
measurement signal shown in Fig. [2] it is instructive to 
plot / and Q as a function of uj m at fixed times t, as 
shown in Fig. [3j With the qubit in red points in 
Fig. O the resonator transmission exhibits the expected 
line shapes for both quadratures. When applying a tt- 
pulse to prepare the qubit in |e), the cavity resonance 
frequency shifts by 2%, but the transmitted quadratures 
respond only on a time scale corresponding to the photon 
lifetime T K = \j k. The lineshape of the cavity transmis- 



sion spectrum centered at -fx wm oru y be reached in the 
limit of T\ ^> T K1 see dotted line in Fig.[3j The interplay 
of the cavity field rise time and the qubit decay time re- 
sults in the observed dynamics of the cavity transmission 
in Figs. [2] and El 

At time t = 180 ns - 1.9 T« ~ 0.2 Ti after the prepa- 
ration of \e) the shift of the cavity resonance to lower fre- 
quency towards -fx is clearly visible, see blue diamonds 
in Fig. El At t = 740 ns - 7.9 T K ~ 0.9 Ti, when « 60% 
of the excited state qubit population P e is decayed, the 
measured curve is approximately the average between the 
steady state \g) and \e) responses, see green crosses in 
Fig. El 

When looking at the time traces in Fig. [2]i), the ef- 
fective shift of the resonance to lower frequency ex- 
plains the reduction of the signal in the Q quadrature 
for measurement detunings A mr > — 0.6 MHz ~ x- For 
A mr < —0.6 MHz the amplitude is increased after the ir- 
pulse because the resonator is driven closer to resonance. 
Given our choice of the rotation of the traces in the IQ- 
plane, the / quadrature of Fig. [2^) acts like a phase and 
always shows a positive response to the 7r-pulse. 

The same considerations explain the features seen in 
the single measurement trace in Figs.^) and c) taken at 
a measurement frequency corresponding to A mr = — x- 
The change of the I and Q quadratures on a timescale 
T K after the 7r-pulse reflects the relaxation of the field 
in response to the qubit excitation. The time scale of 
the return of the quadratures to their initial values is 
determined by the qubit decay at rate 71. 



IV. PULSED MEASUREMENT RESPONSE 

To avoid measurement- induced dephasing during the 
qubit manipulation most of the recent circuit QED exper- 
iments have been performed by probing the qubit state 
with pulsed measurements [27], S, |H, [3J EE El- In 
contrast to a continuous measurement, the measurement 
tone is switched on only after the qubit state prepara- 
tion is completed, see Fig. |4^i) for the pulse scheme. The 
absence of measurement photons during qubit manipula- 
tion also avoids the unwanted ac- Stark shift of the qubit 
transition frequency, thus simplifying qubit control. 

With the qubit in \g) the resonator response reaches 
its steady state at the rate ft, which is seen in the ex- 
ponential rise of the Q quadrature, see blue crosses in 
Fig. HJd). Since the resonator is measured on resonance 
at its pulled frequency A mr = — x, the / quadrature is 
left unchanged, see blue crosses in Fig. St). As in the con- 
tinuous case, the resonator frequency is pulled to uo T + x 
when the qubit is prepared in |e), see red dots in Figs. [3)3 ) 
and c). Since the resonator is now effectively driven off- 
resonantly, the transmitted signal has non vanishing / 
and Q quadrature components both of which contain in- 
formation about the qubit state. With the measurement 
frequency still at A mr = — x, ringing occurs at the differ- 
ence frequency (u r + x) — = 2%. At later times, the 
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Q quadrature [mV] 

FIG. 4: (color online) a) Pulsed measurement scheme. b),c) 
Measurement response versus time for a pulsed averaged mea- 
surement taken in the same conditions as Fig.[2]for the Q and 
/ quadrature respectively. Red points show the response for 
the qubit being prepared in the excited state |e). The blue 
crosses are the measured response to the qubit prepared in 
\g). The trajectory of the field in the /Q-plane is plotted in 
d), the arrows indicate the direction of the time. 



average response is approaching again the steady-state 
value as the qubit decays to \g) at the rate 71. As in the 
continuous case, the qubit lifetime in presence of mea- 
surement photons is obtained from a fit to the Cavity- 
Bloch equations. Note, that the decay of the quadrature 
amplitudes shown in Figs. Ht>) and c) does not directly 
correspond to the exponential decay of the qubit popu- 
lation (<3" z ), but rather is determined by the interplay of 
resonator and qubit evolution. 

The dynamics of the / and Q quadrature amplitudes 
can also be represented in a phase-space plot, see Fig. [Hi). 
The response for the qubit in \g) follows a straight line 
while the response for the qubit in \e) is more complex. 
The nontrivial shape of this curve reinforces the fact 
that both field quadratures contain information about 
the qubit state. It is obvious that a simple rotation in 
the JQ-plane cannot map the signal into a single quadra- 
ture. 

Data taken at different measurement frequencies are 
shown in Fig. [5l As in Sec. the / and Q components 
are rotated such that Q = A and / = in steady-state. 



For the theoretical curves (solid lines), the same set of 
parameters as for the analysis of the continuous measure- 
ment are used, leading to very good agreement. Fig. EK) 
shows the measured response for the qubit in \g). The 
Q quadrature shows the expected exponential rise in the 
cavity population and for t > 0.5 fis we recover the con- 
tinuous measurement response. The I quadrature shows 
the transient part of the response during the initial pop- 
ulation of the resonator, having a negative value (blue 
crosses) for measurements at a frequency above uo v — x 
(blue detuned) and a positive value (red dots) at frequen- 
cies below u Y — x (red detuned). Ringing can be observed 
when the measurement is off-resonant from the pulled 
cavity frequency. Fig. [5b) shows the response with the 
qubit prepared in |e). The response is similar to the one 
shown in Fig. [2] for the continuous measurement, if one 
omits the initial 100 ns where the resonator is populated. 



V. RECONSTRUCTION OF QUBIT STATE 

The detailed understanding of the dynamics of the dis- 
persively coupled qubit /resonator system can be used to 
infer the qubit excited state population p e — ((<r z ) + l)/2. 
Indeed, the difference in the measured response for a 
given unknown state s p (t) and the ground state response 
s g (t), which corresponds to the shaded area indicated in 
Fig. [Dd) and c), is directly proportional to p e . 

To explicitly state this relation, we introduce an ef- 
fective qubit measurement operator M*(t), an approach 
that we have already employed to perform two-qubit 
state tomography using a joint dispersive read-out [28[. 
Here, i = I,Q denote the / and Q field quadratures used 
to measure the qubit state. In terms of this measurement 
operator, the / and Q components of the signal s z p (t) for 
the qubit in state p before the measurement are given by 



<M*(t)> = Tr[pM\t)}, 



(7) 



where M l (t) is determined by the solution to the master 
equation ((4j). Analytical solutions can be found in the 
limit of vanishing qubit decay [Hj], 



M Q (t) 



e Kt ' 2 [2x cos (xt) + ft sin (#)] - 2j( , Q \ 
e — 7^77^ A 8a J 



(k/2)2 



e Kt ' 2 [ttcos(xt) - 2xsin(x*)] - k ^ 



X 2 + (k/2) 2 



which depend on the operator x = A rm + x&z for the 
qubit state-dependent cavity pull. As a consequence 
of performing a quantum non-demolition measurement 
with only a few photons populating the resonator, mix- 
ing transitions between the two qubit states can be ne- 
glected [15| and M l (t) is diagonal at all times. The 
qubit then remains in an eigenstate during the mea- 
surement [37| and we can write M % (t) = s % g (i)\g) (g\ + 

4(t)|e)(e|. The signals s*(t) = Tr[\g) (g\M l (i)] and 

4(t) = Tr[|e)(e|M i (t)] are determined by Eq. ® for the 
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FIG. 5: (color online) / and Q quadratures for pulsed measurements at different detunings from the resonance frequency. The 
plots in a) are taken with the qubit in the ground state \g) while b) displays the response of the system with the qubit prepared 
in the excited state |e). c) shows the result of the pointwise difference of the data acquired with the qubit in the ground state 
and the excited one. The lower panels show time traces taken at different detunings (blue crosses; A mr — 1.4 MHz, black 
diamonds: A mr = 0.3 MHz, red dots: A mr = —0.7 MHz) with comparison to theory (solid lines). 



values (x) = A rm ± x corresponding to the qubit in the 
ground or excited state. To account for qubit relaxation, 
the Cavity-Bloch equations ([5]) are solved to determine 

The qubit excited state population p e (p) in a given 
state p is determined by the normalized area between the 
measured signal s p and theoretical ground state response 

1 silt,) - S*(ti) 

p j p ) = - y — ^4 At, (9) 

where At denotes the discrete time step between data- 
points. s l g / e {tj) are solutions to the Cavity-Bloch equa- 
tions ([5]) with independently determined parameters. Re- 
placing Sp(tj) with the corresponding expressions from 
Eq. (|7J), we notice that Eq. ([9]) simplifies to p e (p) — 
Tr[p|e)(e|], demonstrating that the excited state popu- 
lation of an arbitrary input state is proportional to the 
normalized area between signal and ground state. Thus, 
the effective measurement operator M ,% = \e)(e\ defined 
by this procedure is equivalent to a projective measure- 
ment of the excited qubit state. 

The measurement protocol can be summarized as fol- 
lows: First, the relevant system parameters are deter- 
mined in separate measurements as discussed in Section 
IIIII The qubit lifetime Ti, the single remaining param- 
eter, is determined by applying a 7r-pulse to the qubit 
and analyzing the resulting transmitted signal. From 
this complete set of parameters, the signals s l g (t) and 
s l e (t) are computed. Finally, the excited state popula- 
tion p e is calculated from the recorded signal s l p (t) of an 
arbitrary qubit state p and the theory lines s g (t) and 
using Eq. ©, which amounts to a measurement of 
M' % = |e)(e|. For the particular case of the qubit being 



in |e) after a 7r-pulse, the point-by-point difference signal 
is shown in Fig. (5b). Note, that the excited state popula- 
tion can also be directly inferred from a fit of the Cavity- 
Bloch equations to s p (t) with p e as free fit parameter. 
It is, however, computationally less intensive to calculate 
the population with the area method from Eq. (j9]), that 
is, to perform algebraic operations for the data analysis 
rather than employing a non-linear fit-routine. We have 
checked that both techniques provide the same results 
within the experimental precision. 

To test our method experimentally, we perform a Rabi- 
oscillation experiment [13], where a pulse of variable 
length r and amplitude Q is applied at the effective qubit 
transition frequency cc; S)res . Indeed, the population p e ob- 
tained with the area method (Fig. [6l points) has an rms 




Pulse length [ns] 

FIG. 6: (color online) Rabi oscillations of the qubit popula- 
tion p e reconstructed by analyzing with Eq. © the pulsed Q 
response of the resonator (dots). The black line corresponds 
to the theoretical prediction calculated using Cavity-Bloch 
equations with the following parameters: a; a /27r = 4.504 GHz, 
X /2tt = -1.02 MHz, Ti = 860 ns, Q/2tt = 50 MHz. 
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deviation of less than 1% from the population predicted 
by Eq. ([TDD (Fig. El solid line). The data is also in good 
agreement with a simplified expression 



Pe(t) 



1 1 



-e 4 



|(3 7 i+2 7 0) 



cos(fit/2). (10) 



predicting the time-dependent population of the qubit in 
the limit of large driving fields (ft ^> 71, 7^) [45[. 



in the qubit state is in excellent agreement with calcula- 
tions. The dependence of the measured response on mea- 
surement frequency is well understood both for continu- 
ous and pulsed measurements. From the time dependent 
measurement response we reconstruct the qubit excited 
state population that is used in tomographic measure- 
ments to accurately measure both single and two-qubit 
density matrices [28|, |29[ . 



VI. CONCLUSION 
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